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The oscillator representation is used for the non-perturbative description of vacuum 
particle creation in a strong time-dependent electric field in the framework of scalar QED. 
It is shown that the method can be more effective for the derivation of the quantum 
kinetic equation (KE) in comparison with the Bogoliubov method of time-dependent 
canonical transformations. This KE is used for the investigation of vacuum creation in 
periodical linear and circular polarized electric fields and also in the case of the presence 
of a constant magnetic field, including the back reaction problem. In particular, these 
examples are applied for a model illustration of some features of vacuum creation of 
electron-positron plasma within the planned experiments on the X-ray free electron 
lasers. 
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1. Introduction 

The kinetic description of evolution of particle-antiparticle plasma systems, created 
from vacuum under action of strong fields, has formed now as an active developing 
branch of modern relativistic kinetic theory. Initially, the corresponding Schwinger 
term in the kinetic equation (KE) was introduced on the semi-phenomenological 
basis 21. The first attempts to obtain it on the dynamical level were fulfille d in 
works 121. The fully dynamical ground KE of such kind was obtained in works 
in the framework of QED. Later, large efforts were directed to the application of 
these KEs for the description of pre-equilibrium evolution of quark-gluon plasma, 
arising after collision of two ultrarelativistic heavy ions " , as well as of the forming 
of electron-positron plasma from vacuum in the planned experiments on the X-ray 
free electron lasers (XFEL) at DESY and SLAC ISl. 

However, the existing procedures (Bogoliubov transformation method ^ in com- 
bination with some specific methods of the ki netic theo ry 121 and the formalism of 
the single-time Wigner distribution functions of the dynamical deriva- 
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tion of the KE for the description of the vacuum tunnehng processes are rather 
comphcated and Hmited (especiaUy for particles with many degrees of freedom). 
Originally, the kinetic theory of such kind was constructed on the dynamical basis 
only for the case of charged scalar and Dirac particles (separately and together 
^) in a strong linear polarized space homogeneous time-dependent classical electric 
field (application of these QED equations to the flux tube model in the theory of 
uhrarelativistic heavy ion collision may be based on the Abelian projection method 
EI). On the other hand, large difhculties are met even by the attempt of the transi- 
tion from the electric field of linear polarization to an arbitrary one in the framework 
of the both above mentioned methods. 

Having in mind the transition to the description of more realistic systems, we are 
interested in looking for some alternative, more effective formalism of the derivation 
of the basic KE of the vacuum creation theory. For this aim, we tr y here the oscilla- 
tor (or holomorphic) representation (OR) of quantum field theory l ^^ l ^^ l as a simple 
example of the scalar QED with a strong electric field (Sect. 2). This formalism is 
found very effective in some problem series of the qua ntum f ield theory (e.g., li^). 
The approach based on the Green's functions method l ^^ l ^^ l should be mentioned. 
This method in combination with the Kadanoff-Baym covariant formalism ^1 can 
lead to some alternative variant of the kinetic description of the vacuum particle 
creation. 

In comparison with Bogoliubov method of time-dependent transformations, the 
OR leads at once to the diagonal form of the Hamiltonian (the quasi-particle repre- 
sentation 1^) and to the specific operator equations of motion, which are the basic 
elements for the derivation of the KE El. The KE itself is obtained in Sec. 3 for the 
space homogeneous time-dependent quasi-classical electric field of an arbitrary po- 
larization. As a result of numerical solution of the KE, we show, first of all, that the 
circular polarized harmonic time-dependent electric field is more effective relatively 
to vacuum particle creation in comparison with the case of the linear polarized 
radiation. 

A more complicated case is considered in Sec. 4, where vacuum particle creation 
is investigated in the presence of the time-dependent electric field of linear polar- 
ization and the coUinear constant magnetic field. We study also the back reaction 
problem, based on the union of the KE and the regularized Maxwell equation. Since 
the corresponding generalization of the spinor QED is nontrivial, the examples con- 
sidered in Sec. 3 and 4 can be used to make some qualitative predictions on the 
results of the planned vacuum particle creation on the XFEL experiments. Finally, 
Sec. 5 sums up the results of the work. 

We use the metric g^'' = diag{l, —1, —1, —1) and the natural units h = c = 1 . 

2. The oscillator representation 

We will show as the first step that the holomorphic or the OR l ^^ l ^^ l leads 
at once to the quasi-particle representation in the scalar QED with an external 
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field. So we consider here the case of a complex scalar field in some classical space 
homogeneous time-dependent electric field with 4-potential (in the Hamilton gauge) 



A'^{t) = {0,A\t),A^t),A'it)), (1) 

and the corresponding field strength E = — A , where the overdot denotes the dif- 
ferentiation with respect to time. This field can be considered either as external field 
or as the result of mean field approximation based on substitution of the quantized 
electric field A''{t) with its mean value, {A'^{t)) — A'^{t) , where symbol (...) 
denotes some averaging operation. In the kinetic theory, taking account of fluctua- 
tions leads to collision integrals. Thus, mean field approximation means neglect of 
dissipative effects. 

The Lagrange density 

C{x) = D*^ip*D^'ip - m'^ip*ip (2) 
( Df^ = dfi + ieA^ , e is particle charge with its sign) leads to the equation of motion 



{D^D^ + m^)ip = 0. (3) 

The space homogeneity of the system allows to look for the solution of the Eq. JSJ 
in the following form 

ifix)^ f[dpyi'^ip{p,t), (4) 



where [dp] — (27r) ^/'^d^p . Then the equation of motion of oscillator type follows 
from the Eqs. (O and Q 

(^(±)(p,i)+c.2(p,i)^(±)(p,i) = (5) 
with the time-dependent frequency 



uj{p,t) = \l + P"^ , P = p-eA. (6) 

Symbols (±) correspond to positive and negative frequency solutions of the Eq. Q. 
We suppose a finite limit lim A'^{t) — A'^ in the infinite past and the solutions 

t — * — CO 

ip^^^{p,t) become asymptotically free: 

(/7(±)(p,i) ^ e^''^-*, where w_ = lim uj{t). (7) 

t — ^ — oo t — > — oo 

The Eq. (j^l is the starting point of the definition of time-dependent frequency. Now 
one can introduce the decomposition of the field functions 

/„vpx , ^ 

[dp] { a(-) (p, t) + b^+\-p, t) } , 

V2w(p,i) 1^ J 

/p-ipx , > 
[dp]^== |a(+)(p,t) + b^-\-p,t)] (8) 
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and generalized momenta 



nix) j [dp]^^u;ip,t)e^P^ {a(-)(p,t) - b^+H-P,t)} ■ (9) 

The Eqs. (O and JHl can be obtained from the corresponding decompositions of the 
free field with the help of formal substitution uj{p) = ^ rri^ + p^ u!{p, t) , that is 
admissible for space homogeneous case. 

Thus, the decompositions (jHl) and © are postulated essentially. It is based on 
the possibility of the introduction of the canonical quantization by the complete 
analogy with the case of the free field. Indeed, usual commutation relations for the 
time-dependent creation and annihilation operators 



(-)(p,i),a(+)(g,t)l - \h(-\p,t)M+\q,t)\ ^5{p~q) 



(10) 



etc. follow at once from the canonical commutation relations for the operators <y5(a;) 
and 7r(x) as a direct consequence of the decompositions © and 

Another feature of the OR is connected with the diagonal form of the Hamilto- 
nian in the time-dependent occupation number representation: the substitution of 
the decompositions ^ and © to the full energy of the system 



Ho{t)^ J d^x{Tr*Tr + \/ip*\/ip + m^ip*ip} 



leads to the following Hamiltonian density in the momentum representation 
H{p,t) =uj{p,t) L(+)(p,t)a(-)(p,t)+6(-)(-p,t)fe(+)(-p,t) 



(11) 



(12) 



with u!{p, t) , defined by the Eq. ©. Thus, the OR leads to the quasiparticle repre- 
sentation with the diagonal Hamiltonian density H12|l . This way is found more effec- 
tive in comparison with the " standard" one based on the procedure of Hamiltonian 
diagonalization by the time-dependent Bogoliubov transformation (e.g., ^ '^^ ). 

The OR can be a non-perturbative basis for the solution of one particle problem 
in the strong quasi-classical electric field . So let us obtain the equation of motion 
for the creation and annihilation operators. For this aim, let us also write the action 



J d'^x{Tr*{x)ip{x) +Tr{x)ip*{x) - Ho{x)} (13) 

in this representation 

5 = 1 dtd^pi^^(^a(+Hp,t)d(-Hp,t)-d^+\pM-Hp,t) 
+6(-) (p, <)&(+) (p, t) - (p, i)6(+) (p, t) 
+ A(p, t) (p, i)a(-' i-p, t) - «(+) ip, t)b^+) i-p, t)] ^-Ho (p, t) I , (14) 
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where 

A(P,0 = ^. (15) 

The "anomalous" terms are presented in the Ea. H14|) in the square brackets. These 
terms describe creation and annihilation of the particle-antiparticle pairs with the 
vacuum transition amplitude (|15|l under action of the electric field. Then the op- 
erator equations of motion ^ follow from that after variations by respect to the 
amplitudes a'^^ , 6^^^ and subsequent transition to the occupation number repre- 
sentation with the commutation relations H10|) : 

a(±)(p,t) = iA(p,t)6W(-p,t) + i [i/o(t),a(^'(p,i)] , 

h^^\p,t) = iA(p,t)aW(-p,<) + z [i/o(i),&(±Hp,t)] . (16) 
It is assumed that the electric field is switched off in the infinite past 

lim E{t) = - lim A{t) = 0, (17) 

t — ^ — oo t — ^ — oo 

but the infinitely distant past and future asymptotics of the vector field can be 
distinct i.e. 

lim A{t) = A- 7^ lim A{t) = A+. (18) 

t — > — oo t — >oo 

Let us also notice that the equations of motion (I16II contain some elements of gen- 

Q 99 

eralization in comparison with the known case ■ ' where the field polarization 
remains fixed. This circumstance will be used later on in Sec. 3 in order to de- 
rive the KE for the description of vacuum particle creation in the time-dependent 
electric field of arbitrary polarization. We mark out one feature of the Heisenberg 
equations of motion (|16|l : each of these equations describes some mixture of positive 
and negative energy states ( "non-diagonal" form of equations of motion). That is 
reflected also in the action 1)14(1 . where the "anomalous" contributions ~ a^^^fo'^^^ 
are present. In this sense the OR leads to mixed states. 

The Eqs. H16|l show, that the operator 1)12(1 is not a generator of infinitesimal time 
transformations in the considered Fock space, i.e. it is not the "proper" Hamiltonian 
of the system in this representation. However, the action 1(14(1 is a homogeneous 
square form relative to the operators a^^^ and b^^^ and its anomalous terms can be 
eliminated with the help of some special time-dependent Bogoliubov transformation 

a(-)(p,i) =«t(t)c(-)(p,i) -/3pKt)d'+)(-p,t), 

&(-)(p,i) - a*_p{M-\p,t) - p_pit)cS+\-p,t) (19) 

with the additional relation |ap'(t)p — |/3p*(i)P = 1 (the condition of the transfor- 
mation reversibility and the form-invariance of the canonical commutation relations 
((10(1 \ The substitution of the Eqs. ((1911 into the Eqs. ((16(1 generates to the equations 
of motion for Bogoliubov coefficients : 



ap{t) = i A{p, t)l3t{t)e^^''^P-''^''>\ Pp{t) = 1 A{p, t)at(i)e2^''(^^^*'*°) , (20) 
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where (to is a field switching on time) 



^{p;t,to) = J dt'uj{p,t'). 

to 

Then the action (|14|l is transformed to the standard form 



(21) 



S = / dtd^p{ - 



c(+) (p, t)c<^-^ {p, t) - c(+) (p, i)c(-) (p, t) 



+ S-\p, t)d^+) [p, t) - S-) {p, t)d(+) {p, t)] - H'ip, t) 



(22) 



where the "true" Hamiltonian H'(p,t) has now a non-diagonal form in the mo- 
mentum representation (the exphcit form of H'{p,t) can be found in ^ ). The 
equations of motion 1)16(1 turns into the usual Heisenberg equations, i.e. 



c(±)(p,t) = z[i/'(i),c(±)(p,t) 



(23) 



where the Hamiltonian H'{p,t) is the original point of the theory based on de- 
composition of the type © with the free particle frequency ujo{p) = \/ + p^ . 
Thus, the OR can be considered as the "turned inside out" Bogoliubov method of 
the Hamiltonian diagonalization, 

On the other hand, the Bogoliubov transformation of the type (|19|l can be used 
to restore the diagonal form of the equations of motion 1(161) without mixing the 
terms (the first addends) Then we obtain the following equations of motion in 
momentum representation: 

z^(±)(p,t) = ±f^(p,t)^W(p,t), (24) 

where .4^^^(p, t) are creation and annihilation operators in the new representation 
with some new frequency f2(p, t) . These operators allows us to gain the first integral 
of motion 



{in\A^"-\p,i)A-\p,t)\in) = (m|^(+)(p,0)yl(-)(p,0)|m) 



(25) 



where the averaging procedure is fulfilled over the in- vacuum state. The details con- 
nected with that type representation can be found in the work It is important 
that the correlator 1(25(1 can not be interpreted in quasiparticle terms 



3. Particle creation in the time-dependent electric field 

As an example, we will obtain the KE describing vacuum scalar particle creation 
in the electric field of arbitrary polarization. Let us introduce the distribution 
function of particles using the quasiparticle representation (see Eg. 1(12(1 ') 1^ 

/(p,t) = (m|a(+)(p,t)a(-)(p,i)|m), (26) 

which is related with the distribution function of anti-particles as f{—p, t) = 
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Using the method of works ^ and the basic equations of motion H16|l . it is not 
difficult to get the KE. Differentiating the distribution function H26|l. we obtain 

f{p,t) = ^Aip,t) [/W(p,t) + f^-Hp,t)\ , (27) 
where the anomalous one-particle correlators are 

= {in\a(+\p,t)b^+\-p,t)\tn), 
f^^\p,t) = {in\b'^-\-p,t)a^-\p,t)\in) (28) 
and the vacuum transition amplitude is 

uj^[p,t) 

The equations of motion for the functions (|28|l can be obtained by analogy with 
the Eq. H27|l . We write them out in the integral form 

t 

f^^Hp,t)^^ J dt'A{p,t')[l + 2fip,t')]e^^^'^P-'''''\ (30) 

— OO 

where the asymptotic conditions have been introduced 

lim /W(p,t)-0. (31) 

t — — OO 

The consequence of the Eqs. ijTTI) and ((SU|l is 

f{p,t) = A{p,t)Ref+\p,t). (32) 

Finally, the substitution of the Eqs. H30() into the Eq. I|32() leads to the resulting 
KE 

f{p,t)^I{p,t), (33) 

where the function /(p, t) is the vacuum source term of particles 

t 

Iip,t)^^A{p,t) J dt'Aip,t')[l + 2fip,t')]cos[20ip;t,t')]. (34) 

— oo 

The KE H33|l . H34|l is the non-perturbative result for the electric field in the mean 
field approximation. The particular case of this KE for the linear polarization was 
obtained and researched in detail in works EE1_ 

The KE H33(l can be transformed to a system of ordinary differential equations, 
which is convenient for numerical analysis ^ (to compact our formulas here and 
below, symbols denoting time- and momentum-dependency are dropped in obvious 
cases) 

VI = A(l + 2/) - 2ujV2, 

V2 — 2ujvi, (35) 
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where A is defined by Eq. ^ and vi = iRe/(+' , V2 = ^Im/(+) . The Eqs. §^ 
must be considered with the zero initial conditions. This equation system has the 
first integral 

-il + 2ff + vl + vl = l, (36) 

according to it the phase trajectories are located on two-cavity hyperboloid with 
top coordinates / = wi = W2 = and / = — — V2 = 0. 

Excluding the function / from Eqs. (|35|l . we obtain the two dimensional dy- 
namical system with the equations of motion 

^^i = Ay/l + vl + vl- 2wu2, 

V2 — 2ljv^, (37) 

having non-Hamilton structure. This is the direct consequence of unitary nonequiv- 
alence between the Heisenberg equations of motion H23() (Hamilton dynamics) and 
the Heisenberg type equations l(TB|) . which is the basis of the KE and the 
system H37(l (non-Hamilton dynamics). 

Let us remark also, that the full equation set H35|l have no reversal time sym- 
metry but the physical observables are defined by the functions /, vi only, which 
are not changed at time inversion. 

The system H35(l is integrated via the Runge-Kutta method with the zero initial 
conditions /(p, io) — vi{p,to) = V2{'P,to) = 0. The momentum dependence of 
distribution function is defined by means of digitization of momentum space to 2 or 
3-dimensional grid, in each of its node the system H35|l is solved. The concrete grid 
parameters depend on field strength, the typical values are Ap « O.lm (grid step) 
and Pmax ~ (10 — 15)r7i (grid boundary), so the total number of solved equations 
is of order 10^ - 10^ . 

Now we use the KE ()33|l to compare the effectiveness of particle creation in 
the harmonic time-dependent electric field of circular ( C = 1 ) polarization and of 
linear L£ = 0) one (the analog of last case was considered in the spinor QED in 
works 1^) 

A^(i) = {0, 0, ^x(t) A cos iyt,x{t) A sin vt), (38) 

where x{i) is some function, providing switching on of the field smoothly on a half 
of the field period (the result does not depend on the explicit form of the switching 
function when t ^ 1/^^ )■ Time dependence of particle number density 

nit)^ J[dp]fip,t) (39) 

for the fields ^ with the zero initial conditions at to — is shown in Fig. 1 
(in the natural units). We have chosen the field strength parameter typical for the 
planned experiments on XFEL E, but the considerably larger frequency v — O.lm 
(this is convenient for calculation time reduction). 
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Fig. 1. Time dependence of particle number density for the cases of linear and circular polarized 
field : E = O.lEcr (left panel), E = O.bEcr (right panel). The solid lines correspond to circular 
polarization ( C = 1 ) ^^'^ dotted lines correspond to linear one ( § = ). 



As it follows from Fig. 1, circular polarized field is more effective for vacuum 
particle creation in comparison with linear polarized field of the same amplitude. 
The dynamics of particle creation in periodic field changes qualitatively by the 
increase of the field strength within the range (0.1 — 0.5)Ecr , where Ecr = m? /e . 
Actually, when E = Q.lE^r the time dependence of particle density is periodic 
with the double frequency and approximately constant mean density value for the 
period of laser field (Fig. 1 , left panel) . On the contrary, the mean density increases 
with time when E — 0.5Ecr (the density "accumulation" effect ^3, Fig. 1, right 
panel) . This effect is more expressed in the circle polarized field. We prove that for 
our toy model the creation rate (of charged bosons) in the circle polarized field « 
4 times larger than that of in the linear polarized field. We think that this proper 
effect can be valid for the conditions of the planned experiments on the XFEL. 

If the field strength is of the order of critical one, it is necessary to take into 
account the back reaction of the produced plasma on the primary field. Then the 
external field E{t) should be replaced by the total field 

E{t)^E,n{t)+EeAt), (40) 
where Eex{t) is known function and -Em(t) obeys the Maxwell equation 

% = -e/H5(2/ + .,). (41) 

The total current density in the r.h.s. of this equation is the sum of conductivity 
and vacuum polarization currents. Momentum dependency of the functions /(p, t) 
and vi {p, t) is implicit and it is necessary to investigate their ultraviolet behav- 
ior in order to fulfill the regularization of the integral 1411) . The method of the 
asymptotic expansion can be most adequate one in the framework of the consid- 
ered non-perturbative formalism. The regularization procedure proposed below is a 
variant of n-wave regularization technique Calculation details can be found in 
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the Appendix A. The regularisation in the r.h.s. of the Eq. (|41|l can be achieved by 
substitution vi ^ vi — Vi , where the counter-term Vi{p,t) is the leading term of 
the asymptotic series of the function vi {p, t) over the inverse powers of momentum 

where M is an auxihary mass parameter, M ^ m (som e modi fication of this 
procedure leads to charge renormalization in the Eq. (|41|l ISHEHl). The counter 
terms become significant only in the ultraviolet region. One can choose M 3> A , 
where A is a grid boundary of numerical calculation of the integral H41|) . That 
allows to omit counter-terms as a neglig ible quantity. The KE ^ (or the Eqs. 
H35|l ') and the Maxwell equation H41|) form the complete equation set of the back 
reaction problem, which have been investigated initially in works for the case 
of soliton-like pulse of the external linear polarized field. It was shown that the 
action of the sufficiently strong field E > Ecr excites the plasma oscillations with 
a period, dependent on the residual density after termination of the external pulse. 
In the considered case E < Ecr , the density of excited plasma is small and the 
effects of back reaction are negligible. 



4. Magnetic field influence 

Now we shall investigate the vacuum creation problem in the case of the coUinear 
time-dependent electric and constant magnetic fields 

E(t) = {0,0, E{t)), H={0,0,H). (43) 

That corresponds to the following configuration of vector potential in the Hamilton 
gauge 

A^'{x,t) = {0,-Hxi,0,A{t)). (44) 

Some preliminary results were reported on the conference EH. 

In general case, the introduction of the OR is based on the possibility to define 
the corresponding dispersion law uj{p, t) in the presence of the electromagnetic field. 
This is possible when the spatial and time variables are separable, in particular, 
for the Klein-Gordon equation Q in the field H43|l . That leads to the following 
solutions 

^il.P3(^) - Cr.T^^\t)e-^^/^H^{ii)e'^^^^^+^^^^^ ^ T(±)(t)<I>„,,,,,3 (x), (45) 

where Hn{r]) are the Chebyshev-Hermite polynomials, rj — y^\e\H{xi + p2/\e\H) 
and 

C„ = ^fl|e|Hy^'(2"n!)-i/2 (46) 
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is a normalizing constant. The functions Ti^\t) obey the oscillator-type equation 

fi^'>+Lolit)Ti^'>=0 (47) 

with the dispersion law {P ~ p — eA{t) , p = p^) 

ujl{p,t) ::^m^ + + \e\H{2n+l), ?i = 0,l,... (48) 

The symbols (±) over the functions Tn^^ correspond to positive and negative 
frequency solutions of the Eq. (|17|l at i ^ — oo like in Sect. 3. Finally, the solution 
of of the Eq. H47|l satisfies the relation 

tT<^Ht)dtTi^Ht) ^ Tl, (49) 

as a consequence of the general orthonormalization condition for the solutions of 
the Eq. © 

^ I d^x^<^l^{x) dt [x) = T5nn'5{P2 - P'2)S{P - P'). (50) 

The oscillator decomposition for the field functions and generalized momenta can 
be constructed on the basis of the function set by analogy with the Eqs. © 
and © 

dp2dp 



'^*(a;) = «X! J d-P2dpy^i^nip,t) 

X |a(+)(n,p2,P;i)*n,p2,p(2;) - b^^Hn,P2,P;t)^n,p2A^)} ■ (51) 

The substitution of the Eqs. (|51|l into the action H14|) leads to the equations of 
motion 

d^^Hn,p2,p; t) = ^Wn{p, t)h^^\n, -p2, -p] t) + t[Hn{t), a^^\n,p2,p; t)], 

&W(n,p2,p;i) = ^Wnip,t)a^^Hn,-'P2,-p;t) + t[Ho{t),b^^\n,p2,p;t)], (52) 
where the vacuum transition amplitude is equal 

and Hq is the Hamiltonian in the diagonal quasiparticle representation with the 
energy (|48|l . Let us introduce now the new distribution function of scalar particles 
with momentum {p2,p) on the Landau n -level 

fn{p2,P,t) = {in\a^+\n,p2,p;t)a^^\n,p2,p;t)\in). (54) 
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The resulting KE now has the form 

t t 

fn{p,t)^^Wnip,t) f dt'Wr,ip,t')[l + 2U(p,t')]cOs{2 f dTLU,,{p,T)}. (55) 



2 



It was taken into account here that the distribution function does not depend on 
P2 momentum projection fn{p2,P,t) — fnip,t) due to the absence of the P2 - 
dependence of the amplitude This is the well known effect of the degeneracy 
of the observed quantities on the Landau level We can introduce the effective 
mass TO* = ^y■m^ + \e\H > m in the Eq. l|48|l . which leads to the increase of the 
the effective energy gap width and the suppression of vacuum particle creation. 
This result is in agreement with the kno wn resu lt for a particular case of constant 
electric and magnetic fields combination 1^23 _ 

The KE 155|l can be reduced to the ordinary differential equation system as in 
Sect. 3 

fn = ^WnVin, 
Vin =Wn{l + 2/„) - 2LJ^V2n, (56) 

V2n = 2cj„t;i„. 

The solution of the back reaction problem is based on the KE H55|) and the 
Maxwell equation (the factor \e\H is the consequence of the above mentioned 
degeneracy EHI) 



dE,„ ^ \e\H ^ f , P 
-2e 



n ^ 



fn{P, t) + ^Vin{P,t) 



(57) 



We investigate the back reaction effect at the presence of the critical magnetic 
field H = Her = to^/c = 4.4 • lO^'^Gs for the initial electric pulse of the "rectan- 
gular" form 

Eexit) = (i;,„/2)[tanh(i + b)- tanh(t - b)] (58) 

with the amplitude Em and the duration of 26. It is known that constant 
magnetic field suppress boson creation and enhance fermion creation in the presence 
of the constant electric field. The situation varies essentially for the case of fast time- 
dependent electric field. Figs. 2 and 3 show the evolution of the total electric field 
H40|l and particle number density 

"(*) = 7SE / fniP^t)dp (59) 



{2n) 



n=0 



for the weak and strong electric field with H = Her and = 10 , here e is the 
effective colour charge (charge value selection is motivated by the flux-tube mod el 
in the theory of the QGP generation at ultra-relativistic heavy ion collision EEl^. 
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Fig. 2. Back reaction in the weak electric and strong magnetic field { H = Her )'■ total electric 
field (left panel) , particle number density (right panel) . The stable plasma oscillations arc excited 
after termination of the external pulse when E > Ecr ■ 



Total field almost coincides with external pulse with = O.lEcr , i-e. the back 
reaction is negligible here due to the small particle density. The internal field be- 
comes appreciable when — 0.5Ecr and stable plasma oscillations are excited 
after the termination of the external pulse. The period of oscillations decreases with 
the growth of the residual particle density as ^/n . The amplitude of the oscillations 
becomes large-scale modulated at the further increase of the field, Fig. 3. 

During the action of the external pulse (|58|l with the flat top, particle density 
grows linearly. The rate value (without the magnetic field) is well approximated 
with the Schwinger formula for the field strength close to critical value but in the 
weak field these quantities differ strongly. Analogous situation is for the suppression 
coefficient of boson creation which is equal to sinhx/x for E = const, where 
X = ttH/E I^. For the top plot on Fig. 2 this value is « 2 • 10^^ , whereas the real 
difference obtained in the kinetic approach is about 50 % . 
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t*m t*m 

Fig. 3. Back reaction in the strong electric and magnetic field; total electric field (left panel), 
particle number density (right panel). The plasma oscillations becomes large-scale modulated at 
the increase of the field to E K. E^r • 

5. Summary 

We have shown on an example of scalar QED, the oscillator (holomorphic) repre- 
sentation is now, apparently, the most effective non-perturbative method for the 
kinetic description of vacuum particle creation. For an illustration, we have ob- 
tained the KE for the case of the arbitrary time-dependent space homogeneous 
electric field and applied it for a comparison of the effectiveness and other features 
of vacuum tunneling processes in the linear and circular polarized fields. 

Since the corresponding calculations in the spinor QED were not yet fulfilled, 
these investigations can be the foundation for the qualitative understanding of pro- 
cesses of electron-positron plasma creation in strong laser fields in the planned 
DESY and SLAG installations. We have investigated also the influence of the mag- 
netic field on vacuum particle creation and have studied the corresponding back- 
reaction problem, where we have observed great deviations with respect to predic- 
tion of the well known results for the case of the constant electric and magnetic 
fields combination 02]. 

The results obtained here can be interesting also for the development of the 
fiux-tube model of superconductive type, describing the pre-equilibrium evolution 



1, 2008 16:42 WSPC/INSTRUCTION FILE Osc6ra 



The Kinetic Description of Vacuum Particle Creation... 15 



of quark-glupn plasma, generated under conditions of ultrarelativistic heavy ion 
collisions The present approach reveals some perspective in the investigation 
of the dynamics of vacuum particle creation in the strong space homogeneous back- 
ground fields in more realistic situations when one deals with vacuum creation of 
particles with inner degrees of freedom (electron-positron, quark-gluon plasma etc.). 
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Appendix A. Regularization procedure 

The ultraviolet behavior of the unknown functions / , vi and V2 is determined by 
the dispersion law uj{p,t) and the vacuum transition amplitude A(p, i) . We can 
write them out in the spherical coordinate system and define their asymptotic in 
the ultraviolet region for the electric field ^ 

eE{t)nq{t) eE{t)n 1 

Ap,i = ^ 7T{ 777: (A.l 

uj^{q) q^oo q{t) q{t) 



where q{t) = y P^(t) and n — P/q . Now we construct the asymptotic series over 
the inverse powers of momentum q 

/ = Vs = J2vi-\ s = l,2, (A.2) 

n— 1 n—1 

where /("^ , vi^^ ^ q~" . We get the following leading contributions using these 
decompositions and the Eqs. (|35|l 

rU)f ,N [ eE{t)n \ ^ (3) eE(t)n (2) eE{t)n 

(9'^) ^ 1 ) ' "1 ="1 (9'^) = «2 (9,0 = ^^- (A.3) 

The unique purpose of these counter terms is the elimination of ultraviolet 
divergence in the integrals defining the densities of physical quantities (e.g., Vi 
provides the regularization of the electric current in the r.h.s. of the Eq. (|41|l . - 
of the energy density etc.). That is why it is necessary to introduce some mechanism 
of switching on of the counter-terms ljA.31) only in the asymptotic region q ^ 00 . 
Switching on the counter terms in the region q ~f 00 can be achieved, e.g., by the 
replacement of q"^ ^ q"^ + in the Eqs. (A3), where is an auxiliary mass 
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parameter (Pauli-Willars procedure), and trending il/ — > oo after calculations of 
the integral. Another possibihty is to modify the counter terms, e.g., 

f^£>{q,t) = f^\q,t)FM{q,t) (A.4) 

with the switching function of the type FM{q, t) = I — cxp(— g/Af) . 

When using the computer calculations and therefore introducing the cut-off 
momentum parameter (A = qmax is a grid boundary), one can always assume 
counter terms to be negligible small and omit them. 
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